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ABSTRACT 


This  thesis  investigates  the  problem  of  aerial  search  for  a  mobile 
surface  target.  The  initial  position  of  the  target  is  known,  and  a 
certain  amount  of  time  elapses  before  the  search  begins.  The  target 
has  a  variety  of  choices  as  to  speed  and  direction,  and  the  search  air¬ 
craft  has  a  choice  of  search  patterns  to  be  flown.  The  method  proposed 
to  solve  the  problem  of  maximizing  or  minimizing  the  probability  of 
detection  is  the  game  theoretic  analysis  of  various  strategies  avail¬ 
able  to  the  searcher  and  the  evader.  Four  strategies  for  each  are 
compared  and  evaluated,  with  optimal  strategies,  within  the  set  of 
strategies  investigated,  being  derived. 
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I.  Introduction. 

A  frequently  encountered  problem  of  attack  aircraft  is  the  detection 
and  destruction  of  mobile  enemy  installations  at  some  time  after  initial 
detection.  For  example,  a  flight  may  be  fired  upon  by  a  mobile  surface- 
to-air  missile  site.  Upon  return  to  the  air  base  or  carrier  the  location 
of  the  attack  is  reported  and  a  strike  against  it  is  readied  and  sent  out. 
In  the  meantime,  the  missile  site  has  often  had  time  to  depart  the  posi¬ 
tion  in  which  it  was  originally  situated.  There  are  many  other  ways  in 
which  a  mobile  target  may  be  detected,  such  as  aerial  reconaissance,  air¬ 
borne  electronic  detection,  or  visual  sighting.  Some  of  these  may  give 
a  probabilistic  rather  than  a  deterministic  position. 

The  mobile  target,  if  it  knows  it  has  been  discovered,  has  the  op¬ 
tions  of  remaining  in  its  present  position  or  moving.  If  it  moves,  it 
has  some  variety  in  its  choices  of  direction  and  speed.  Often  these 
choices  of  movement  are  to  some  degree  limited.  The  search  aircraft  are 
then  confronted  with  the  problem  of  finding  and  attacking  the  moving 
strike  objective.  This  can  be  hampered  by  such  variables  as  time  to  sun¬ 
set,  terrain,  weather,  and  other  uncontrollable  elements.  The  target  may 
be  found  by  methods  other  than  visual  if  it  is  emitting  electronic  sig¬ 
nals,  i.e.,  if  it  has  not  moved  or  if  it  has  moved  and  is  set  up  and 
operating  in  another  position.  The  most  difficult  problem  arises  when 
the  method  of  detection  must  be  strictly  visual,  as  when  the  objective 
is  moving  or  if  it  has  set  up  in  another  position  and  is  not  emitting. 

The  degree  of  difficulty  also  depends  on  how  much  freedom  of  motion  the 
target  has.  If  the  direction  of  travel  is  limited  rather  strictly,  as 
when  there  are  only  a  small  number  of  roads  along  which  it  may  travel 
and  off-road  travel  is  not  possible,  the  searcher's  task  is  somewhat 
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simplified.  If,  however,  it  is  free  to  travel  in  any  or  almost  any  direc¬ 
tion,  the  problem  becomes  much  more  difficult.  This  thesis  examines  some 

’  i'  *  •  .  ; 

of  the  problems  confronting  a  searcher  faced  with  the  problem  of  discover- 
ing  a  target  with  unrestricted  direction  of  movement,  and  with  varying 
maximum  speeds. 

This  problem  may  be  attacked  in  two  fashions.  One  method  is  to  con¬ 
sider  the  movements  of  the  target  as  being  determined  probabilistically, 
with  such  variables  as  direction,  velocity,  and  time  of  movement  all  being 
described  by  some  density  function,  so  that  an  estimate  of  its  most  likely 
position,  or  most  likely  set  of  positions,  can  be  found  using  the  usual 
rules  of  probability.  In  fact,  all  these  variables  have  some  probabilis¬ 
tic  limitations,  but  within  these  limitations  there  is  another  influence: 
tr.e  judgment  of  the  commander  of  the  site.  It  would  be  somevrtiat  naive  to 
think  that  the  commander  would  be  guided  solely  by  the  constraints  imposed 
by  the  state  of  nature.  Rather,  it  is  more  reasonable  to  assume  that  he 
would  try  to  vary  his  behavior  within  the  limits  of  his  physical  con¬ 
straints  so  as  to  attain  maximum  liklihood  of  evasion.  If  this  did  not 
occur  immediately,  then  certainly  experience  would  lead  him  to  adopt  more 
sophisticated  evasive  tactics.  The  intelligent  searcher  would  then  alter  s 

his  search  pattern  to  counter  the  evasive  tactics  to  as  great  a  degree  as 
possible.  This  development  brings  the  problem  into  the  area  of  game 
theory,  and  it  is  with  this  aspect  that  this  paper  concerns  itself . 

To  simplify  the  problem  a  number  of  probabilistic  variables  are 
treated  deterministically.  The  effect  of  weather,  time  of  day,  terrain 
and  variance  in  terrain  over  the  search  area,  and  other  external  influ¬ 
ences  are  ignored.  These  and  other  assumptions  are  discussed  in  the  next 
section.  A  small  number  of  strategies  for  both  searcher  and  evader  are 
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discussed.  While  an  attempt  has  been  made  to  consider  the  most  advan¬ 
tageous  strategies,  there  are  probably  many  of  merit  which  have  not  been 
included.  It  is  felt  that  an  approach  has  been  made  to  the  solution  of 
the  problem,  however,  and  some  possibly  beneficial  tactics  derived. 

Hie  general  approach  is  to  select  several  tactics,  both  for  the 
searcher  and  the  target.  These  tactics  are  then  evaluated  in  pairwise 
combinations,  using  game  theory  criteria.  This  allows  conclusions  to  be 
drawn  concerning  the  relative  merits  of  the  strategies. 

The  next  section  presents  the  basic  assumptions  of  the  model  and 
attempts  to  justify  them.  Hie  ensuing  sections  develop  the  general 
theoretical  framework  and  the  specific  mathematical  formulations  of  the 
proposed  strategies.  After  comparing  these  strategies,  some  conclusions 
as  to  the  value  of  each  are  inferred. 


II.  Theoretical  Basis  for  the  Model. 

2.1  Assumptions. 

It  would  be  an  almost  impossible  task  to  obtain  a  universal  solution 
to  this  search  problem,  in  which  all  relevant  factors  were  included.  On 
the  other  hand,  elimination  of  a  critical  variable  would  leave  the  con¬ 
clusions  worthless.  The  assumptions  made  in  the  development  of  the  model 
are  set  forth  here  in  an  attempt  to  explain  its  limitations  and  clarify 
its  objectives. 

The  space  in  which  the  game  is  played  is  a  two-dimensional  Euclidean 
geographical  space.  The  terrain  problem  is  ignored  in  that  the  area  is 

« 

considered  to  be  perfectly  flat,  and  in  no  way  affects  either  target  mo¬ 
bility  or  the  searcher's  visibility.  Hie  target  is  free  to  move  in  any 
direction  and  at  any  velocity  up  to  its  maximum  for  any  duration  of  time. 
The  searcher  has  a  fixed  range  sighting  probability;  within  this  range  the 
probability  of  detection  is  1.0,  outside  this  range  the  probability  of 
detection  is  0.  The  searcher  is  capable  of  flying  any  pattern  exactly. 

The  initial  position  of  the  target  is  known  exactly.  The  target  knows 
he  has  been  detected,  and  takes  action  accordingly.  No  variables  such 
as  weather,  time  of  day,  or  other  external  phenomena  are  considered.  The 
time  at  which  the  aircraft  will  arrive  to  begin  his  search  is  known  to 
the  target  as  soon  as  the  initial  detection  occurs. 

There  are  other  minor  assumptions  in  the  model.  These  cited  are 
considered  to  be  the  most  relevant-  To  be  able  to  claim  that  the  model 
has  any  relevance  to  the  real  world  some  of  these  assumptions  must  be 
justif led. 

First,  consider  the  playing  area.  Unless  the  terrain  is  extremely 
rough,  the  aircraft  acts  essentially  as  if  it  were  flying  above  a  two- 


dimensional  area.  The  behavior  of  the  target  is  more  affected  by  this  as- 
sunption.  Since  the  target-chosen  factors  are  speed  and  direction,  the 
effect  of  terrain  on  these  is  relevant.  In  a  real-world  situation,  the 
target  might  indeed  be  somewhat  terrain  limited.  However,  most  missile 
and  radar  sites  are  near  cities,  which  in  turn  are  usually  situated  in 
flat  terrain.  Also,  there  is  usually  a  heavy  density  of  roads  near  any 
city,  which  gives  the  target  a  great  degree  of  flexibility  in  its  move¬ 
ments,  both  as  to  direction  and  speed.  If  the  area  around  the  initial 
position  is  not  travelable  in  almost  any  direction,  the  problem  is  simpli¬ 
fied  to  some  extent,  but  is  susceptible  to  the  same  type  analysis  as  is 
proposed  here. 

The  assunption  that  the  target  can  move  at  any  speed  up  to  its  max¬ 
imum  for  the  duration  of  the  search  is  reasonable  if  the  search  aircraft 
arrive  within  a  reasonable  length  of  time,  on  the  order  of  four  hours 
after  initial  detection,  and  stay  on  station  a  moderate  length  of  time. 
Since  most  attack  circraft  can  fly  at  low  altitudes  for  only  about  two 
hours,  the  total  time  the  target  must  be  able  to  keep  moving  is  about 
six  or  seven  hours.  This  is  well  within  the  capability  of  most  mobile 
sites. 

The  method  of  first  detection  of  the  target  is  often  visual  or  photo 
graphic,  in  which  cases  the  position  is  known  exactly.  If  the  position 
is  not  known  exactly,  the  point  of  highest  probability  is  usually  taken 
as  the  true  position  and  treated  deterministically  in  any  actual  search. 
The  effect  of  large  errors  in  the  position  of  the  original  site  is  not 
examined  in  this  discussion. 

The  implications  of  a  fixed  range  of  detection  are  well  known  and 
no  further  implications  have  arisen  in  this  analysis.  [ 3 ]  However,  the 
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idea  that  an  aircraft  can  fly  any  desired  pattern  exactly  will  strike  any 
experienced  aviator  as  a  flagrant  misconception.  In  fact,  this  is  prob¬ 
ably  the  most  critical  assvinption  of  the  model.  Recent  years  have  seen 
the  development  of  Doppler  navigational  equipment  which  enables  an  air¬ 
craft  to  fly  a  precise  track,  and  airborne  radar  in  a  control  aircraft 
can  guide  other  aircraft  over  a  desired  pattern  with  good  accuracy.  There 
is  no  reason  to  believe  that  such  techniques  could  not  be  applied  to  search 
aircraft  if  the  situation  warranted. 

Another  critical  assumption  is  that  the  target  knows,  as  soon  as  de¬ 
tected,  at  what  time  the  search  aircraft  will  arrive.  Some  examination 
of  the  actual  procedures  reveals  this  assumption  to  be  less  restrictive 
than  it  initially  appears.  If  the  target,  moves,  the  probability  of  de¬ 
tection  diminishes  with  time  in  any  of  the  strategies  considered.  The 
search  aircraft  will  therefore  want  to  begin  the  search  as  rapidly  as 
possible.  The  target  often  knows  the  position  of  the  base  or  carrier 
from  which  the  attack  aircraft  will  be  launched,  and  can  therefore  pre¬ 
dict  fairly  accurately  their  arrival  time.  An  estimate  might  also  be  ob¬ 
tained  from  observing  past  arrival  times. 

While  these  assumptions  restrict  the  application  of  the  model  to 
a  limited  situation,  they  allow  the  formulation  of  a  general  mathematical 
model.  Refinements  of  the  general  approach  may  be  made  to  deal  with 
specific  deviations. 

2.2  Outline  of  Game  Theory  Approach. 

The  searcher,  denoted  player  A,  and  target,  denoted  player  B,  have 
available  to  them  a  large  and  possibly  infinite  set  of  strategies.  Each 
attempts  to  maximize  his  payoff  by  choosing  judiciously  among  these 
strategies.  B  may  choose  any  direction  and  any  speed  up  to  his  maximum, 
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and  may  stop  his  movement  at  any  time.  A  may  choose  any  search  pattern 
over  the  playing  area  which  does  not  exceed  his  range.  The  payoff  func¬ 
tion  for  A  is  the  probability  of  target  detection.  The  payoff  to  B  is 
either  the  negative  of  the  probability  of  detection,  -P(D),  or  the  prob¬ 
ability  that  A  does  not  detect  B,  l-P(D).  If  the  former  is  defined  as 
B's  payoff  the  game  is  zero-sum;  if  the  latter  the  game  is  constant-sum 
and  the  sum  is  one. 

The  set  of  strategies  for  player  A  is  denoted  f1  and  each  element  of 
this  strategy  set  is  denoted  •<  .  The  set  of  strategies  for  player  B  is 
denoted  ^  and  the  elements  of  the  set  are  denoted  fb  .  The  pairs  of 
strategies  («S/3)  are  elements  of  the  Cartesian  product  of  P  and 
denoted  H  x  ^  .  The  payoff  to  A  for  a  given  pair  of  strategies,  <7S> 
and  ,  is  M(^0  ,  fi0  )  =  P(A  detects  B/  o<  •  n  *  fi 0  ).  For  an  opti¬ 
mal  strategy  to  exist,  there  must  exist  strategies  e**  and  fi*  such  that 

min  max  M(®<  ,  fi )  *  max  min  M(o <  ,/S)  -  M  (®<  , /3  )  A7 

P  <*  °<  fi 

Whether  in  fact  such  strategies  exist  has  not  been  determined.  To 
determine  the  existence  of  such  strategies,  either  all  possible  strate¬ 
gies  must  be  examined,  or  it  must  be  proven  topologically  that  such  strat 
egies  exist.  The  first  method  is  difficult  due  to  the  large  number  of 
possible  strategies.  There  exist  some  existence  theorems  which  pertain 
to  a  pursuit-and-evasion  game,  but  none  were  found  which  were  applicable 
to  the  game  under  discussion.  An  existence  theorem  developed  by  Ryll- 
Nardzewski  [l]  which  was  the  most  general  as  to  existence  criteria,  indi¬ 
cated  that  a  solution  involving  pure  strategies  did  not  exist. 

It  is  also  possible  that  an  optimal  solution  which  is  a  convex  com¬ 
bination  of  various  strategies  exists.  To  find  this,  it  is  necessary  to 
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examine  the  payoff  matrijc  of  the  various  strategies,  and  to  apply  ttie 
usual  methods  of  game  theory, to  find  a  solution.  This  requires  develop¬ 
ing  the  probability  of  detection  for  various  strategies.  Therefore  it  is 
worthwhile  to  develop  payoff  functions  for  individually  considered  strate¬ 
gies  of  A  and  B. 

Tne  general  procedure  then  is  to  choose  some  strategy  oi0  and  some 
strategy  .  fl>e  specifies  a  rule  for  choosing  direction,  speed  and 
stopping  point  for  the  target.  <y#  specifies  a  pattern  to  be  flown  over 
the  playing  area.  The  payoff  is  P(A  detects  b  /o<.  ^  fl  ~  )•  These 

strategies  may  be  probabilistic  and  specify  a  probability  function  for 
choosing  speed  or  initial  aircraft  search  point,  or  may  be  deterministic, 
such  as  always  start  the  search  at  a  radius  of  five  miles,  or  always  re¬ 
main  at  the  initial  position.  If  the  set  of  attractive  strategies  is 
compared  pairwise  in  this  manner,  the  matrix  of  strategy  payoffs  is  gen¬ 
erated  . 
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III.  Development  of  Strategies. 

3.1  Strategies  Considered. 

As  has  been  noted,  there  are  a  great  many  possible  strategies  for 
both  searcher  and  target.  After  considering  a  large  number  of  strategies, 
the  following  were  selected  for  investigation  on  the  basis  of  their  intui¬ 
tive  attractiveness  or  because  they  are  frequently  used  in  practice.  A 
glossary  of  terms  is  given  in  table  one.  The  strategies  as  proposed  by 
the  author  are: 

■V 

Player  A  searches  over  v  in  the  most  efficient  manner  which  also 
searches  over  all  directions  of  possible  target  travel,  (0,  2fr-).  This  is 
taken  to  be  a  spiral  search,  and  the  choice  of  the  target  velocities  to 
be  searched  is  determined  in  a  game  theoretic  manner. 

Table  I.  Glossary  of  Terms 

vfc  target  speed 

vfl  search  aircraft  speed 

vt  maximum  target  speed 
max 

Q  direction  of  target  movement  or  polar  angular  coordinate 
of  aircraft  position 

R  aircraft  range 

T  .  minimum  time  after  sighting  at  which  target  may  move 
min 

T  time  from  initial  detection  of  aircraft  arrival 
arr 

tn  T  -  T  .  or  time  the  target  has  been  in  motion  at  aircraft 
0  arr  min 

arrival  time 

r  radius  from  initial  target  position 

rQ  radius  at  which  the  aircraft  starts  its  search 
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circumference  of  the  circuit  searched  by  the  aircraft 
W  fixed  sighting  range 

A  diagram  of  a  spiral  search  is  given  in  Figure  1. 

C<  . 

2' 

This  is  essentially  the  same  as  OC  with  the  additional  condi¬ 
tion  that  the  original  target  position  is  always  investigated. 

C<  . 

3 

The  aircraft  search  randomly  in  the  circle'of  radius  (v.  *  t), 

max 

or  the  circle  within  which  the  target  must  be. 

c<  . 

4* 

The  aircraft  search  uniformly  over  v  and  (0,  2 /T  ) 

The  evasive  strategies  are 


The  target  chooses  velocity  and  direction  to  counter o<  ^  accord¬ 
ing  to  a  game  theory  criteria. 
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3.3  Spiral  Search. 

It  is  possible  to  demonstrate  that  a  spiral  search  for  a  target 
moving  radially  outward  from  a  fixed  point  is  the  most  efficient  way 
of  searching  over  various  target  velocities,  if  the  target  direction  is 
uniformly  distributed  on  (0,  2jh  ) .  This  may  be  intuitively  justified  as 
follows:  Suppose  a  searcher  has  a  given  fixed  range.  He  desires  a 
search  pattern  which  will  minimize  his  path  length  while  not  searching 
over  any  velocity  range  more  than  once.  For  a  single  circuit,  i.e.,  a 
search  of  9  from  0  to  2 ft,  the  target  velocities  depend  upon  search  width 
W  and  the  initial  radius  of  the  search,  r^.  If  the  velocity  range  is  to 
remain  constant  throughout  the  circuit,  the  search  radius  must  move  out¬ 
ward  with  the  velocity  which  would  have  put  the  target  at  r^,  which  is 
r0/t0.  If  the  searcher  moves  outward  in  such  a  manner,  and  increases  the 
radial  component  of  his  velocity  with  each  subsequent  circuit  such  that 
rR  m  rn  ^  -  (r^  /  tQ)  *  t,  he  will  not  re-search  any  previously  searched 
velocity.  This  pattern  may  be  flown  in  several  ways.  The  outward  radial 
component  of  velocity  may  be  chosen  to  remain  uniform  throughout  the 
search.  Another  method  is  to  vary  the  radial  component  after  the  comple¬ 
tion  of  each  circuit,  thereby  flying  a  series  of  spirals  joined  by  short 
cross-legs.  While  this  is  slightly  less  efficient  than  a  continuous 
spiral,  it  has  the  advantage  that  it  searches  over  a  continuous  target 
velocity  range.  This  makes  it  more  compatible  with  an  analysis  of  the 
optimum  velocity  choices  of  the  target,  and  for  this  reason  this  method 
is  used. 

The  equation  for  a  single  circuit  is  derived  as  follows: 

The  search  aircraft  flies  an  increment  of  the  circuit  of  lengthA  S 
in  timeA  t.  During  this  time  the  target  moves  radially  outward  with 
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The  final  forms  of  the  expressions  are 


S 


1  ) 


r  =  rQ  e 


v10 

a 


0 

v 

To  get  these  expressions  in  a  more  manageable  form,  express  e  a 


as  an  infinite  series. 


He 

va 


2 

4  02 

1  +  Vt  8  +  v? 
v. 


a 


+  .  .  .  • 


a 
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For  vt  »  vfl, 


lie 

e  va  i 


1  +  — —  0 
v 

a 


Since  vfl  is  usually  over  300  and  is  usually  less  than  30,  this 
is  a  valid  approximation. 

Then 

va  vt 

S  =  — _  r  (  i  +  — —  0  -  1  )  =  3r 
vt  °  va  ° 


and 

vt 

r  =  r  (  1  +  — L-  e  ) 
va 

At  the  end  of  each  circuit,  the  searcher  must  position  himself  to 
search  over  a  new  velocity  increment.  If  the  detection  range  either  side 
of  the  flight  path  is  W,  the  minimum  velocity  unsearched  is  (rQ  +  W)/  tQ, 
which  is  that  velocity  which  would  put  the  target  at  Tq  +  W  at  the  start 
of  the  search.  Since  the  search  width  is  2W,  a  radius  of  rQ  +  2W  will 
search  velocities  from  (rg  +  W)  /  tg  to  (r^  +  3W)  /  tQ,  if  the  velocity 
is  assumed  to  be  constant  across  the  search  width.  This  is  not  strictly 
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r0  +  W  r0  +  2w  r0  +  3  u 

true,  since  - - -  <  -  . .  -  <  —  .  .  ;  however  since  this  dis- 

to  to  to 

tance  is  small  it  will  be  ignored.  For  example,  if  W  =  1,  tQ  -•  2,  then 
WAq  is  .5,  which  is  the  maximum  deviation  from  the  velocity  at  the  cen¬ 
ter  of  the  search  strip.  As  the  time  increases,  this  deviation  becomes 
smaller. 

If  the  target  were  at  Tq  +  2W  at  tg,  it  would  move  radially  outward 
during  the  time  the  first  circuit  was  being  flown.  The  time  required  to 
search  the  first  velocity  increment  is 


Then  the  initial  radius  for  the  following  circuit  is 


rn  +  2w  Sn 

rl  =  r„  .  2W  .  -  (-2.) 

*0  va 


The  distance  from  the  position  at  the  end  of  the  first  circuit  to  the 
beginning  point  of  the  next  circuit  is 


)  =  2W  + 


4W  r. 


2  r0 

r  *  ri  “  r0  (  1  +  -^-r - 

al0 


If  r_  =  100,  W  =  1,  v  =  300,  and  tn  =  2, 


r  =  2  +  8„  'A0n°  -  12.16 

2’  300 

Since  these  cross-legs  are  small  compared  to  total  aircraft  range,  they 
will  be  ignored  in  further  computations. 

Then  to  compute  any  r,  it  is  necessary  to  take  its  position  at  tg 
and  to  add  to  it  the  distance  traveled  while  the  aircraft  was  searching 


previous  circuits. 


rn  =  r0  +  2nW  +  (- 


rn  +  2nW 
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since 


n-1 


n-1 

Because  of  the  term  ^  S,,  a  recursive  relation  exists. 

i=0 


A  solution  for  N  *  is  possible  if  the  total  path  flown  is  re¬ 

stricted  by  the  aircraft  range,  or 


N 


i=0 


However,  because  of  the  complexity  of  the  expression  an  analytic  solution 

is  very  difficult  and  an  iterative  method  of  solution  for  the  maximum 

velocity  searched  as  a  function  oj.’  rQ  is  laborious.  For  these  reasons  a 

computer  program  was  used  to  solve  the  expression.  Thus  a  method  exists 

for  solving  for  the  range  of  target  velocities  searched  for  any  rQ,  with 

input  parameters  W,  R,  v  ,  vt  and  tn.  This  solution  is  necessarily 

a  max  u 

discrete,  since  the  rg's  are  discrete;  however,  any  degree  of  accuracy 
desired  is  obtainable. 

3.4  Game  Theory  Evaluations  of  the  Proposed  Strategies. 

3.4.1  Discussion  of  Strategyo<  . 

Using  the  previously  derived  expression  for  the  range  of  velocities 

searched  over  for  any  r^,  the  minimum  and  maximum  target  velocities  may 

be  obtained.  The  probability  of  detection  for  any  given  initial  search 
★ 

radius,  say  rQ,  is 

P  [v  .  (r*)  <  v  4  v  (r*)  /  r  ■  r  *  / 

i  min  o  t  max  o'  I  o  ol 


Then 

M(eV  Z3 i>  *  p  (°) 

a  P  ( v  •  (r*)  <  v  <  v  (r*)  |  r  =  r*  J  P  (t  -  r*  ) 

i  min  o  t  max  of  o  o/  \o 

The  lower  the  value  of  r*  the  greater  is  the  range  of  target  vel¬ 
ocities  searched.  For  the  last  circuit,  generally  not  all  values  of  0 
are  searched  since  the  circuit  must  terminate  when  maximum  range  is  reached. 
The  target  direction  of  travel  is  distributed  uniformaly  on  (0,  2  ),  which 

means  the  probability  of  detection  on  the  last  circuit  is  the  percentage 
of  the  circumference  searched.  This  probability  is  denoted  9^.  The  den¬ 
sity  function  for  P(D/rQ  *  Tq)  is 


The  graph  of  this  density  function  is  depicted  in  Figure  2. 

Denote  the  range  of  velocities  searched  by  vr-  The  payoff  function 

is  M  (v*,  v*)  =  P(  D  /  v  =  v*  and  v..  =  v*). 
r  t  r  r  t  t 

The  distribution  functions  for  v_  and  v*  are  F(v  )  and  G(v.).  These 

r  t  r  t 

are  not  fixed  distributions,  but  are  any  possible  distributions.  The  ex¬ 
pected  value  of  the  game  for  two  distributions  F  and  G  is  denoted  E(F,G). 

vt  vt 

rmax  /max 

E(F,G)  =  J  J  M(vr,  vt)  dF(vr)  dG(vt> 

0  0 
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Figure  2 

Conditional  Density  Function  of  the  Probability  of  Detection 

If  M(v  ,  v  )  is  continuous  then  a  solution  is  guaranteed.  However,  in 
r  t 

the  case  of  spiral  search  M  is  not  continuous,  and  therefore  no  solution 
is  guaranteed.  However,  an  examination  of  the  proof  of  this  theorem  in 
4  indicates  that  the  requirement  of  a  continuous  payoff  function  may  be 
too  restrictive.  The  assumption  of  continuity  is  necessary  in  the  proof 
to  guarantee  the  existance  of  the  Stieltjes  integral  of 


f M(x,  y)  dF(x) 
M(x,  y)  dG(y) 
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and  to  guarantee  the  existance  of  certain  minima  and  maxima  over  the  do¬ 
main  of  the  payoff  function.  It  has  not  been  determined  if  the  existence 
of  a  solution  to  the  problem  under  discussion  exists. 

The  existence  of  a  game  theory  solution  has  been  presumed  and  den¬ 
sity  functions  for  the  searcher  and  target  derived.  Due  to  the  complex¬ 
ity  of  the  analytic  expression  giving  the  range  of  target  velocities 
searched  for  a  given  r^,  an  approximation  to  the  solution  using  discrete 
values  of  Tq  and  vt  has  been  used. 

If  Tq  and  vfc  are  allowed  to  take  on  discrete  values  and  the  probabil¬ 
ity  of  detection  derived  for  each  combination  of  these  discrete  values,  a 
matrix  approximating  the  continuous-domained  payoff  function  may  be  ob¬ 
tained.  The  degree  of  accuracy  is  dependent  on  the  fineness  of  the  incre¬ 
ments  of  v^  and  rQ;  for  purposes  of  the  present  investigation  increments 
of  one  unit  of  speed,  i.e.,  one  knot  or  one  mile  per  hour  were  used.  The 

resultant  matrix  is  m  x  m+1,  where  m  =  v.  .  There  are  several  meth- 

max 

ods  to  obtain  the  game  theory  solution  of  this  matrix,  but  the  precise 
methods  of  solution  have  the  detriment  of  being  extremely  complex  for  a 
matrix  of  any  size.  For  this  reason  an  algorithm  which  gives  an  approx¬ 
imate  solution  was  used.  The  algorithm  and  the  computer  program  used  to 
perform  the  calculations  are  discussed  in  Appendix  I.  An  extensive  dis¬ 
cussion  with  proof  of  convergence  is  given  in  Karlin  C2J.  The  relevant 
properties  of  the  method  are: 

(1)  An  approximate  strategy  for  A  and  B  may  be  found. 

(2)  Upper  and  lower  bounds  on  the  value  of  the  game  are  generated. 

The  accuracy  of  the  solution  depends  on  the  nunber  of  iterations  per¬ 
formed.  Two  thousand  iterations  are  considered  to  give  a  usable  approx¬ 
imation  for  the  matrix  under  discussion.  As  an  accuracy  check  one  program 
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was  run  which  performed  twenty  thousand  iterations.  The  upper  and  lower 
values  of  the  game,  denoted,  V  and  V,  were  .1254  and  .0845  for  2000  iter¬ 
ations  and  .1153  and  .1017  for  20,000  iterations.  The  total  computer 
time  time  for  the  more  accurate  method  is  four  minutes  and  seven  seconds. 
With  2,000  iterations  ten  sets  of  input  parameters  may  be  run  in  three 
minutes  and  fifty  seconds.  Using  this  program  optimal  strategies  for  A 
and  B  were  found,  with  A  restricted  to  using  only  spiral  search.  The 

strategies  for  searcher  and  target  for  varying  v.  and  t  are  tabulated 

tmax  0 

in  Appendix  II. 

3.4.2.  Discussion  of  Strategy®^- 

All  of  the  theory. and  general  discussion  relating  tooi  ^  is  appli¬ 
cable  toO<  2*  The  same  computer  program  with  appropriate  modifications 
was  used  to  solve  the  payoff  matrix.  The  range  of  the  search  aircraft 
was  diminished  by  the  radius  of  initial  search  to  compensate  for  search¬ 
ing  vt  0  on  each  search  pattern,  which  had  little  effect  on  the  out¬ 
come.  This  strategy  forced  B  to  always  move  from  the  position  of  initial 
detection,  which  to  some  extent  alleviates  the  assunption  that  the  time 
the  target  begins  its  movement  is  a  deterministic  rather  than  a  probabil¬ 
istic  quantity.  The  reason  for  this  assumption  is  that  there  is  a  lower 
bound  on  the  length  of  time  the  site  takes  to  disassemble  components  and 
makes  ready  to  move.  If  B  is  forced  to  move,  the  longer  he  is  in  motion 
the  more  difficult  is  the  problem  of  detection.  If  B  knows  the  initial 
position  will  always  be  searched,  he  is  more  likely  to  make  an  effort  to 
become  highly  efficient  in  preparing  to  move,  and  therefore  the  range  of 
times  at  which  he  moves  cluster  more  closely  around  the  minimum  move  time. 
The  solution  generated  by  the  computer  for  this  strategy  indicates  that  B 
will  move  at  a  very  low  speed  in  about  the  same  proportion  that  he  did  not 
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move  when  A  plays  strategy 

3.4.3.  Discussion  of  Strategy®^ 

This  strategy  requires  A  to  search  randomly  in  a  circle  of  radius 

(t0  ■  v.  ) .  It  has  the  advantage  of  not  requiring  any  expenditure  of 
max 

effort  on  navigation,  and  in  fact  is  often  used  in  practice.  It  is  an 
inefficient  search  method  and  ia  included  largely  as  a  basis  of  compari¬ 
son  for  more  effective  strategies.  The  area  in  tdiich  it  searches  is  not 
constant  if  the  target  is  moving  radially  outward,  and  the  probability 
of  detection  is  less  than  is  given  in  the  following  formula.  However, 
since  the  probability  of  detection  is  low  for  this  method,  no  accuracy 
was  lost  in  the  general  conclusions. 

The  probability  of  detection  is  given  by 

-2WR 

PCD)  =  1  -  e  A  £3] 

The  derivation  of  this  formula  is  given  by  Koopman.  [ 3 ] 

3.4.4.  Discussion  of  Strategy0*^. 

In  this  strategy  A  searches  all  velocities  of  B  with  equal  likeli¬ 
hood,  or  vf  is  distributed  uniformaly  on  (0,  v  ).  There  is  a  problem 

max 

in  this  strategy  arising  from  the  fact  that  searching  some  velocities 
requires  less  distance  traveled  than  others.  The  strategy  is  therefore 
formulated  as  the  sequential  choice  of  velocities  to  be  searched  with 
each  choice  being  uniformly  distributed  on  the  target  velocity  range. 

The  expected  probability  of  detection  based  on  the  expected  range  of 
velocities  searched  is  used  in  the  payoff  matrix  of  the  various  strate¬ 
gies. 

3.4.  5.  Discussion  of  Strategies  and  ^8 

These  strategies  are  defined  in  the  course  of  deriving  strategies 
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and  o<  2’  As  with  ^  and&<2,  discrete  values  are  used  to  approximate 
the  continuous  strategy. 

3.4.6.  Discussion  of  Strategy  y 

B  moves  at  time  tn  with  a  velocity  chosen  uniformly  on  (0,  vt  ) 
u  max 

and  continues  to  move  with  this  velocity  throughout  the  time  period.  The 

direction  of  travel  is  chosen  uniformly  on  (0,  2W).  This  distribution 

is  also  obtained  if  the  target  moves  with  velocity  v  and  stops  at 

rmax 

times  uniformly  distributed  on  (0,  tg)  where  0  denotes  the  minimum  time 
at  which  the  target  may  move. 

3.4.7.  Discussion  of  Strategy^  . 

In  this  strategy  player  B  distributes  himself  uniformly  over  the  geo¬ 
graphical  area  of  radius  v  .  tg.  The  probability  of  detection  is 

max 
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and  the  distribution  is 

P(vt  <  vr) 


Since  strategies  ^  and  have  been  approximated  using  discrete  vel- 
oc.ity  increments  this  strategy  may  be  defined  in  velocity  increments  using 


the  distribution  function,  which  gives  an  expression  compatible  with  •< 
ando<  2- 


IV.  Evaluation  of  Strategies. 

4.1.  Method  of  Evaluation  of  Strategies 

The  pairwise  comparison  of  all  combinations  of  strategies  considered 
may  be  made  and  the  results  tabulated  in  a  matrix  of  payoffs.  Examina¬ 
tion  of  this  matrix  then  reveals  the  optimum  or  combination  of  optimum 
strategies.  The  elements  of  the  matrix  are 


*ij  ■  pfD  **  •  "*  i  °  P  =  0 J 

The  general  expression  for  the  probability  of  detection  for  strategies 
at  and  o<  2 


i=l  i=l 


max  uu 

P(D)  =  £  P(r0  =  r £)  P(D  /  rQ  ■  r  ) 


^P(v  .  (r.)  v  v  (r.)) 

V  mm  J-  t  max  1 

The  comparison  of  ^  and  ^ ,  and  of  Qs  ^  and^  2  ^ave  been  made  in  the 
solution  of  the  game.  The  solution  for  elements  a12  and  a21  ma;  be  made 
using  the  general  formula  directly.  To  solve  for  elements  a1  ,  a  ,  a  , 
824  the  following  modification  of  the  basic  formula  is  used: 

vtmax  Vmax  (ri> 

P(D)  *  £  P(r0  »  ri)  J  t(vt)  for  fi  ^  and  ^ 

1-1  vmin  <ri> 

where 


/  - 


vz  -  . 

max  mm 


for  /5 


;ma 
f(vt)  =  — 


vmax  “  vmin 


for 


29 


for  a  single  circuit. 


2TT  h 


max 


Therefore  a^  was  taken  to  be 


(E(C)) 


The  payoff  matrix  for  parameters  R  ■  500,  tQ 
and  V2  *  300  is 


2,  W  *  1,  v  «  30 
Lmax 


1 

2 

3 

4 

.108 

.103 

.180 

.071 

.103 

.101 

.174 

.074 

.041 

.041 

.041 

.041 

.068 

.064 

.077 

.054 

This  matrix  has  a  saddle  point  at  a2^.  Therefore  pure  strategies  exist 
for  players  A  and  B.  It  appears  that  searching  by  strategies  oC  1  and 
o<2  highly  desirable  for  player  A,  since  the  values  of  ajj  and  a2j 
are  very  close  and  a  deviation  by  B  from  strategy  fi  ^  greatly  increases 
the  probability  of  detection.  If  B  distributes  himself  uniformly  over 
his  range  of  velocities  then  A  increases  the  probability  of  detection 
by  more  than  a  factor  of  two  if  strategies  •<  j  and  OC  2  are  used.  To 
determine  tfiether  the  existence  of  a  saddle  point  in  this  case  was  unique, 
values  of  the  elements  of  the  first  two  rows  of  the  matrix  were  calcu¬ 
lated  for  six  other  sets  of  input  parameters  and  in  each  case  a  saddle 
point  existed  at  either  a^  or  a2^. 
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For  0(  ^  the  formula  for  the  probability  of  detection  of  a  uniformly  dis¬ 
tributed  target  is  used  although  this  is  not  true  for  strategies  ft 

an<*  ft  g*  However,  the  low  probability  of  detection  for  random  search 
indicates  that  more  precise  computation  is  not  necessary. 

For  fi  4  there  is  a  difficulty  in  that  path  length  for  searching  dif¬ 
ferent  velocities  is  not  constant.  This  has  been  handled  by  using  the 
expected  value  of  the  path  length  for  a  single  circuit.  The  expected 
value  of  v  is  calculated  for  each  strategy  of  B  and  this  is  used  to  get 
the  expected  nunber  of  circuits.  This  is  then  used  to  calculate  the  prob¬ 
ability  of  detection. 

The  expected  target  velocity  is 

vt 

Lma/ 

E(v  )  =  21  v*  p(v4.  =  v  ) 

0 


Then 


Thus 


E  £path  length]  =  e[L]  =  2  E(vt> 


E  ^number  of  circuits ]  *  e£c]  *  — j-. 


While  it  would  be  possible  to  solve  for  the  exact  probability  of  detec¬ 
tion  of  0<  4  for  fll  i  and  2  this  was  not  done  since  rough  calculations 
indicated  that  such  a  value  would  be  very  low.  The  expected  nunber  of 
circuits  may  be  used  to  solve  for  a^j  directly  by 

2W 

E(C)  x  tQ 


max 


since  the  target  speed  is  distributed  uniformaly  on  (0,  vt  ).  The  ex- 

lmax 

pected  probability  of  detection  for  a^  has  been  calculated  to  be 
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V.  Conclusions. 


5.1.  Evaluations  of  Strategies 

The  low  value  of  strategies  ^  ®nd  the  random  and  uniform  on 

vt  searches,  indicate  that  their  use  in  an  actual  situation  would  not  be 
optimal.  However,  the  use  of  a  simple  strategy  that  corrects  some  of  the 
deficiencies  of  these  strategies  may  be  used.  If  an  area  search  is  used 
which  does  not  re-search  any  area,  the  probability  of  detection  will  be 
approximately 
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This  is  very  slightly  less  than  the  probability  of  detection  using  strat¬ 
egies  OC  and  the  sP*-ral  searches,  and  is  very  much  simpler  to  apply. 

It  must  be  concluded  that  a  uniform  geographic  distribution  of  the 
target  is  optimal  in  the  set  of  strategies  discussed.  Since  no  counter 
to  this  strategy  seems  to  exist,  this  may  well  be  an  optimal  solution 
over  all  strategies  available  to  player  B.  If  the  equation  for  probabil¬ 
ity  of  detection  used  for  non-overlapping  flight  path  over  a  circle  of 

radius  v  *  tQ  is  used,  a  very  close  approximation  to  the  probability 
max 

of  detection  given  by  the  saddle  point  is  obtained.  This  indicates  that 
strategies  o*  ^  and  &<  2  are  optimal  in  the  set  of  strategies  studied  by 
virtue  of  their  characteristic  of  non-overlapping  flight  paths  rather 
than  because  they  derive  an  advantage  because  of  some  other  characteris¬ 
tic.  If  this  is  true  and  if  B  plays  strategy^}  uniform  geographical 
distribution,  the  effort  needed  to  apply  strategies  ^  and  o<  2  seems 
unjustified.  However,  the  distinction  between  uniform  geographical  dis¬ 
tribution  and  uniform  velocity  distribution  is  not  obvious,  and  conceiv¬ 
ably  a  naive  player  B  could  choose  to  play  strategy  /&  ^  and  distribute 
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uniformly  on  (0,  tt  )•  Since  using  strategies  *  ,  or  <*,  in  this  case 

max 

«uld  considerably  increase  the  probability  of  detecting  the  target,  the 
anticipated  benefits  might  outweigh  the  increased  effort.  This  might  be 
augmented  by  an  attempt  to  estimate  the  type  of  strategy  applied  by  the 
evader  by  goodness-of-f it  tests,  to  determine  the  probability  that  B  rs 
using  some  strategy  other  than/3  „  and  to  use  a  counter-strategy  which 
takes  maximum  advantage  of  B's  mistakes. 
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APPENDIX  I 


Computer  Program  for  Solving  for  Strategies  1  and  2 
This  program  has  two  basic  parts.  The  first  develops  the  range  of 
target  velocities  searched  for  various  values  of  Tq  using  a  spiral  search, 
and  generates  a  payoff  matrix  using  increments  of  one  unit  of  speed  be¬ 
tween  successive  values  of  Tq  and  vt.  The  second  part  of  the  program 
solves  the  matrix  using  a  method  for  approximating  the  solution  of  matrix 
games . 

The  inputs  are: 


RO 

VA 

TO 

RNG 

DELR 

W 

TMAX 

The  outputs 
VMIN 
VMAX 
THETA 
RO 

A 

X(i) 

YCi) 


initial  radius  searched 
search  aircraft  speed 

length  of  time  the  target  has  been  in  motion  when 
search  begins 

search  aircraft  range 

amount  the  initial  radius  is  incremented  for  each 
successive  search 

search  width 

maximum  target  speed 

are: 

minimun  target  speed  searched  for  a  given  RO 

maximum  target  speed  searched  for  a  given  RO 

percent  of  circle  searched  on  last  circuit 

radius  corresponding  to  VMIN,  VMAX  and  THETA  on  same 
line  of  printout 

matrix  of  detection  probabilities 

approximation  of  optimum  percent  of  time  the  search 
should  start  at  RO  equal  to  i  •  tg  +  W 

optimum  percent  of  time  the  target  should  choose 
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VI  upper  value  of  game 

V2  lower  value  of  game 

Internal  symbols  ares 

SUM  totpl  distance  flown 

S  c ire  inference  of  present  circuit 

SLAST  percent  of  circumference  searched  on  last  circuit 

VT(j)  target  velocity  equal  to  (j~l) 

vt 

'-max 

SUMA  (IA)  ^  j  for  j  .)  SUMB(j)  is  maximum  over  j 

i=l 


SUMB(IB) 


.  ^  SUMA(i)  is  minimum  over  i 


NA  minimum  SUMB  (j) 

NB  maximum  SUMA  (i) 

YB  (IE)  number  of  times  the  target  has  picked  velocity  YB 
(IE  -  1) 

XA  (IX)  number  of  times  the  search  aircraft  has  picked  initial 
velocity  searched  XA®  IX  up  to  present  iteration. 

This  program  is  written  in  Fortran  63  for  use  on  the  CDC  1604  com¬ 
puter.  The  first  part  of  the  program  uses  the  equation  developed  in 
Section  3.3  to  generate  minimun  and  maximum  target  velocities  searched 
for  various  values  of  R0.  The  initial  R0  is  an  input  to  the  program,  as 
is  the  increment  of  R0,  designated  DEIit.  For  purposes  of  this  investi¬ 
gation  increments  of  one  unit  of  speed  were  desired,  and  therefore  the 
statement  DELR  =  TO  has  been  inserted.  R0  normally  begins  at  W,  which 
assures  search  of  the  initial  target  position.  Each  successive  R0  de¬ 
fines  a  row  of  the  payoff  matrix.  After  VMIN  and  VMAX  are  found,  the 
spectrum  of  target  velocities  are  subjected  to  IF  tests  to  find  if  they 
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fall  within  or  without  the  detected  velocities.  They  are  assigned  prob¬ 
abilities  of  detection  of  0  or  1.  Then  the  element  of  the  row  correspon¬ 
ding  to  VMAX  is  assigned  a  probability  of  detection  THETA.  This  process 
is  continued  until  RO  =  VTMAX  *  TO.  At  this  point  the  matrix  is  com¬ 
plete. 

Part  II  of  the  program  obtains  an  approximate  game  theory  solution 
for  t  is  matrix.  Obtaining  an  exact  solution  for  a  large  matrix  is  diffi¬ 
cult  unless  a  saddle  point  exists,  which  requires  that 

max  min  a—  =  min  max  a., 
i  j  j  i 

But  for  each  j  in  the  matrix  there  exists  an  element  a. .  such  that 
J  1J 

a^..-  =  0.  For  each  i  there  exists  an  a.  .  such  that  a.  .  0.  Therefore 

no  saddle  point  exists. 

The  exact  methods  of  computation  involve  examination  of  all  adjoints 
of  the  submatrices  of  the  matrix.  The  generation  of  the  adjoints  re¬ 
quires  evaluation  of  the  inverse  and  the  determinant  of  the  submatrix. 

The  number  of  submatrices  of  a  30  x  30  matrix  is  There- 

n=l 

fore  a  method  approximating  the  actual  solution  has  been  used.  This  uses 
an  iterative  process  in  which  player  A  chooses  his  strategy  by  looking 
at  the  column  vector  which  is  the  sum  of  all  columns  chosen  by  B,  and 
playing  the  row  corresponding  to  the  maximum  element  in  this  vector.  B 
looks  at  a  row  vector  composed  of  the  sum  of  all  rows  chosen  by  A  and 
chooses  the  minimun  element  for  his  next  play.  The  number  of  times  each 
strategy  is  chosen  divided  by  the  total  number  is  an  approximation  of 
the  optimal  strategies.  The  minimum  element  in  the  row  vector  and  the 
maximum  element  in  the  column  vector  divided  by  the  number  of  iterations 
give  lower  and  upper  bounds  on  the  value  of  the  game. 
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The  row  vector  of  sums  is  SUMB  in  the  program,  and  the  column  vector 
of  sums  is  SUMA.  The  number  of  times  each  row  and  column  has  been  chosen 
are  XA  and  YB.  Statements  70  through  80  choose  the  maximum  over  SUMA. 
Statements  74  through  77  choose  the  minimun  over  SUMB. 

The  program  as  presented  generates  the  solution  for  strategy*#^. 

To  obtain  the  matrix  and  solution  for  strategy j,  statement  97  and  the 
preceding  statement  may  be  deleted,  along  with  the  fourth  statement  after 
statement  17. 
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